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^ ' Abstract. All covariant time operators with normalized probability distribution are 

derived. Symmetry criteria are invoked to arrive at a unique expression for a given 
Hamiltonian. As an application, a well known result for the arrival time distribution 
of a free particle is generalized and extended. Interestingly, the resulting arrival 
time distribution operator is connected to a particular, positive, quantization of the 
£T) \ classical current. For particles in a potential we also introduce and study the notion 

C"*** ' of conditional arrival-time distribution. 
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1. Introduction 

X. 

In spite of a famous footnote by Pauli JI] on the purported non-existence of a time 
operator in quantum mechanics, physically motivated questions relating to time occur 
quite naturally in the form of arrival times or of time durations. These become 
interesting when the extension of the wavefunction comes into play, and a quantum 
mechanical description by operators is not at all obvious. Also clock time operators have 
been investigated in the literature. For reviews of several aspects of these difficulties 
and of some more recent results see, e.g. [21 E]. As quite recently stressed by us and 
Munoz in [I] , for a given system there is no single time operator. Indeed, there may be 
many, corresponding to different observables and measuring devices. 

There are two main groups of physically relevant time operators. These are 
associated with time durations and time instants, respectively. An example of a duration 
is the dwell time of a particle in a region of space. The corresponding operator commutes 
with the Hamiltonian since the duration does not depend on a particular reference 
instant |5j. In the other group, the time observables are covariant in the sense that 
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when the preparation time is shifted they are shifted by the same amount, either forward 
(quantum clocks) [6] or backward (event times recorded with a stopwatch, e.g. time of 
arrival), and they are conjugate to the Hamiltonian. 

The measurement of a clock time operator T is supposed to yield, up to a constant, 
the parametric time. More precisely, T should satisfy (ipt\T\ipt) = t+ (ipo\T\ipo) . Time-of 
arrival operators are mathematically similar to clock time operators, but the physics is 
different. Some of their properties are motivated by the behaviour of a classical particle, 
and in any case may be put at test experimentally. A free particle in one dimension 
must, during its motion, arrive at a given point with certainty, and it must also arrive 
in three dimensions at an infinitely extended plane. On a half-line, a free particle is 
reflected at the end point and may therefore arrive twice at another point. In this 
case it is reasonable to investigate first (or second) arrivals at a point |4J. The most 
difficult case is a particle interacting with an external potential [7J EJ [9j QUI EED, EE] . The 
particle can be deflected and, in one dimension, reflected so that it may, with a non-zero 
probability, never arrive at a particular position. In such a case it is then natural to ask 
for the probability distribution of first arrival at time t, given that the particle is arriving 
at all, i.e. to ask for the conditional first-arrival probability distribution, which is, by 
definition, normalized to 1. In jl] the most general form of covariant and normalized 
time observables was announced and a proof of a special case was indicated. The general 
proof of this result will be given in this paper. Reference [I] also gave applications to 
arrival times for a free particle on a half line and to Lyapunov operators in quantum 
mechanics. Based on physically motivated postulates pertaining to arrivals at a plane 
in three dimensions, [13] and [H] derived a time-of-arrival operator for a free particle 
in three dimensions coming in from one side of the plane. Among the postulates were 
invariance of the time operator under Galilei transformations, which leave invariant the 
plane in question, and minimal means-square deviation. This result will be generalized 
in the present paper and the underlying conditions will be weakened. In particular, 
Galilei invariance will not be needed. 

The plan of the paper is as follows. In Section |2] we study covariant and normalized 
time observables and explicitly spell out their most general form, with the detailed 
proof referred to the Appendix. In Section [3] we analyze physical conditions, in 
particular symmetry conditions, under which time observables become unique, and 
discuss applications to clock-time operators. In Section @] we give applications to time- 
of-arrival operators for free particles in one and three dimensions. In particular, the 
result of [T3~| [H] will be generalized. The restriction on the incoming direction will be 
lifted and the above Galilei invariance will be replaced by invariance under translations 
which leave the plane in question invariant. It will be shown that the resulting arrival 
time distribution and operator are related to a particular, positive, quantization of the 
classical current. Due to the backflow effect [15] the usual quantization in the form of 
the quantum mechanical probability does not have the necessary positivity properties. 
In Section |5] we study the notion of first arrivals in the presence of a potential and the 
associated probability distribution operator, which in general is not normalized to 1. 
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Normalizing by the total arrival probability for a given state one obtains the conditional 
arrival probability distribution which is normalized, but, by construction, is not bilinear 
in the wavefunction. Alternatively, to keep bilinearity we consider normalization by 
means of an operator and the resulting distribution operators. A one-dimensional 
example for a particle in a potential symmetric about the origin and arrivals at the origin 
is discussed. In the Appendix mathematical fine points such as regularity properties are 
discussed. 



2. Clock time and time-of-arrival operators 

Clock time operators can be associated with a quantum clock whose average measures 
the progressing parametric time. A good clock would also minimize the variance in order 
to estimate the time as accurately as possible with a finite number of measurements. 

For a clock observable, let the probability of finding the measured time in the 
interval (oo,r) for a state \ip) be given by the expectation of an operator F T , which 
plays the role of a cumulative probability operator, so that F^ = 1. The derivative 
d{ip\F T \ip) / dr is the corresponding temporal probability distribution and it is normalized 
to 1. For a momentum measurement the analogous operator would be a projector, while 
here one just has < F T < 1 with selfadjoint F T . The collection of F T 's yields a positive 
operator-valued measure. We define the probability distribution operator fl r and the 
time operator T by 

n T = —f t , 

ar 

f = J drrfl T . (1) 

For details on the derivative see the Appendix. The mean value of observed time can 
be written as 

<V| J dTril T \ii) = ^\T\^) . (2) 

The second moment, if it exists, is given by (0| J drT 2 fl T and similarly for higher 
moments. 

Covariance of a clock time operator with respect to ordinary (parametric) time 
means that the probabilities of finding the measured time in the interval (— oo, r) for 
state |^o ) an d in (— oo,r + t) for \tp t ) coincide. This implies 

F T = e~ iflT/h F e iflT/R . (3) 

From this one obtains by differentiation (see the Appendix for details) 

n = y[#,F ], (4) 
fL T = e- iAT/h n e iAT/fi . (5) 

Note that flo and fl T are in general not operators on the Hilbert space but only bilinear 
forms evaluated between normalizable vectors from the domain of H. Since {ip\F T \ijj) 
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is non- decreasing, fl T is a positive bilinear form. An expression like (E\fl Q \E f ) has to 
be understood as a distribution. Since the diagonal E = E' has measure it is no 
contradiction that (jlj) formally gives on the diagonal while explicit examples may give 
a non-zero value [4]. By a change of variable in ([1]) one obtains 

e iHt/hf e -iHt/h = f + t (q) 

Furthermore, H and T satisfy the canonical commutation relations as bilinear forms on 
a suitable domain. 

A superscript A denotes quantities for (first) arrival times at a given location. In 
contrast to clock times, if the particle's state is shifted in time by to, it should arrive 
a time to earlier, and the temporal probability distribution should be shifted by to to 
earlier times. Thus the analog of the cumulative probability operator in ([3]) and the 
probability density operator must now satisfy 

pA = e im/hpA e -iAt/R (7) 

n£ = i (8) 

f t A = e i6t/hf s A e -iAt/h^ (9) 

Again, il^ is a positive bilinear form. 

In case of a particle in a potential, it may happen that for a given state the particle 
arrives at the location with a total probability less than 1. This means that the integral 
of Ilf is not unity. For the average arrival time one can then not use the left-hand side 
of the analog of ([2]) but rather 

drrfl T \7p)/(^\N\^), (10) 

where N = J drfl T = F^. However, if the total arrival probability is always 1, one can 
define an arrival time operator by 



T A = I drrU T (11) 
dtte iflt ' h f^e-' lflt / h . (12) 



As seen by a change of integration variable, an arrival time operator behaves as the 
negative of a clock time operator whose probability distribution operator is given by 
fl T = fl^ T , with a cumulative probability operator F T = 1 — F^ T . 

We now prove the following general result for covariant clock operators and arrival 
time operators with normalized distributions. This was already announced in [4] without 



proof. In Appendix A.l and Appendix A. 2 we present rigorous mathematical details. 



Let if be a Hamiltonian with purely continuous eigenvalues and eigenvectors 
\E,a), where the degeneracy index can be assumed to be a discrete number. Their 
normalization is taken as 

(E, a\E', a'} = 5 aa <5(E — E'). (13) 
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For simplicity we assume the same degeneracy for each E. We consider any covariant 
clock time operator which has a probability distribution normalized to 1 and whose 
second moment exists for a dense set of vectors. Then for t = the probability 
distribution II is given by 



2^5^/ d E dE'^2h{E,a)\E,a){E',a'\bi{E',a') 



(14) 



where bi(E, a) = (E, oc\ 6j) are any functions satisfying the two conditions 

^2b i (E,a)b i (E,a')=5 aal , (15) 

i 

\dEbi(E,a)\ 2 integrable over finite intervals. (16) 

i 

The same results hold for Hq of normalized arrival time distributions. The complete 
distributions fl t and flf and the associated time operators are then given by ([!])- (f5ll 
and (jTj) - ffTTj) . respectively. 

For a state with ^(i?, a) = (E, a\ip), one has 



(17) 



dE^i)(E,a) ihd E ^{E,a) 

a 

+ J dEj2^( E , a ) a ') bi ( E > a "> ihd Eh(E, a') 

aa' i 

second moment = h 2 J dE \ c)e h(E, a) ip(E, o 



A mathematically detailed proof will be given in |Appendix A. 1| and |Appendix A. 2 



For simplicity we indicate here a formal proof for the case that the (continuous) 
eigenvalues of H are non-degenerate. We start by constructing the functions 6j. For a 
given n , one can choose a maximal set {|<7i)} of vectors satisfying 

(ftlftotoi) = • (19) 

Such a maximal set is easily constructed by the standard Schmidt orthogonalization 
procedure. Then a possible set of |£>j)'s is given by 

|&i> = V27rHn o |0i) . (20) 
Indeed, from ( 1201) one has 

^-^2(get\k)(bi\g p ) = ^2(g a \flo\g i )(g l \flo\g^ 

i i 

= Sgjdpi = 8 a p. (21) 



Equation (TT9j) then implies 

no = ^E^)^l' ( 22 ) 
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which is (Q3J). It should be noted that the |6j)'s in the decomposition of IIo m (|14p are 
not unique. 

The \bi) 's need not be normalizable and the bi(E) 's could in principle be 



distributions. That the bi(E) 's are in fact functions will be shown in Appendix A.l 



They have to satisfy certain properties in order that the total probability is 1 and that 
the second moment is finite. Indeed, for given normalized state the total temporal 
probability is, with ip(E) = (E\ijj), 

— (V|e- iff */ ft n e iH */ fi ^) 
= E/^/ dEdE'e-^ E - E '^ h ^)b t (E)bW)^E') 
= I dEdE'S(E - E')i^b i (E)k{E 1 )'iP(E') 

= £/ dEj{E)Y,HE)HE)^E) = l . (23) 

Since \ip) is arbitrary this implies 

J2k(E)k(E) = l. (24) 



The formal use of the 5 function will be justified in Appendix A.l In a similar way one 
obtains 

(i(,\f\ij))= J dE$(E)ihip'{E) 

+ j dE\^{E)\ 2 Y,HE)ihbW) . (25) 

Note that Yl^iK * s purely imaginary, from (1241) . and thus vanishes if 6, is real. The 
second moment is 

— t 2 (ij\e- iHt / h fl e [ ^ h \ij) 
2n 

H J2j^J dEdE'dEdve-^-^^b^bW)^') 
h 2 Y,jdE d E (vp) k(E) ) d E (k(E)iP(Ef 



= n 2 dE{\i//w\*+Y,\Km 2 \i>m 3 

+ 2ReY^HE)y i (E)iJ{E)i;'(E)} (26) 

% 

by (j24p . This is finite if the contribution from the first and second term are finite, 
and for the latter to hold for all infinitely differentiable functions ip(E) with compact 
support one must have that YliWi(E)\ 2 is integrable over any finite interval. Again, 
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differentiability of the fej's will be shown in the Appendix, as well as the proof in the 
case that H has degenerate eigenvalues. 

Conversely, it is easily checked that any functions bi(E,a) which satisfy ( fl5l) and 
( Tl6|) give rise to a covariant time operator with normalized probability distribution and 
second moment for a dense set of vectors. 



3. Degenerate energy eigenstates: unique clock time operators from 
symmetry and minimal variance 

Clock time operators have a normalized probability distribution, and this will be 
assumed in the following. To reduce the multitude of covariant clock time operators 
we may require that the variance AT, or mean-square deviation AT 2 , is minimal for all 
states for which the second moment exists. Physically this means that no other time 
observable can be measured more precisely. However, requiring minimal variance by 
itself does not make T unique, not even in the case of non-degenerate spectrum of H, 
as was shown in Ref. [3]. But if in addition one restricts the set of functions 6j by 
symmetry requirements, uniqueness may then follow from minimality of AT, as will be 
illustrated now. 

We denote the anti-unitary time reversal operator by 9. In the x space 
representation one has (0?/>)(x) = ?/>( x )- If the Hamiltonian is time reversal invariant, 
we may demand 

ef(S) = -f (27) 

and similarly for the probability distribution. By ([3]) this implies 

efl o = n o . (28) 

In [3J it was shown for the non-degenerate eigenvalue case that time reversal invariance 
of the Hamiltonian H and minimal AT together imply that T and fl t are unique and 
given by 

ft* = dEdE'e-^- E 'W h \Ee)(E' e \, 

f = J dttflt, (29) 

with time reversal invariant \Eq). 

If the Hamiltonian has degenerate continuous eigenvalues this is no longer true and 
one needs additional conditions to obtain uniqueness, as discussed now. 

Example 1: In one dimension we consider 

H = P 2 /2m + V(\X\) (30) 

and assume a purely continuous spectrum. H is invariant under space reflection, denoted 
by a, and under time reversal 6. We assume T and its probability distribution to have 
analogous properties. For ft this implies 

<3-n <7 = en e = eafi d-e = n . (31) 
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Then, under the assumption of normalizability, the time operator becomes unique, as 
shown now. As eigenstates of H we choose real and either symmetric or antisymmetric 
wavef unctions, denoted by \E, ±) and normalized as 5{E — E'). Then Q\E, ±) = \E, ±) 
and a\E, ±) = ±\E, ±). Decomposing the functions bi a in (fi~4")) for n into real and 
imaginary parts one finds from flo = |(fto + ©fto©) that (lT4"j) can be written with real 
functions when one uses \E, +) and \E, — ) as basis. Writing again il = f (Ao + crflocr) 
one sees that the mixed aa' terms cancel so that one can write Ho in the form 

ft °=2^EE I dEdE'b ia (E)\E,a)(E',a\b ia (E') (32) 

i a=± J 

with real functions b ia . Note that (I3"2l has no mixed aa' terms and is in this sense 
diagonal. The normalization condition (fl5|) is replaced by 

J2b ia (E) 2 = l , a = ±. (33) 

i 

Diagonality in a and reality of the functions b, ia will be crucial for determining b ia 
explicitly. In fact, expanding the right-hand side of (fl8|) for the second moment, the 
term 



(34) 



J dEb' i (Ea)b i (E,a') ip(E, a) ip'(E, a') + c.c. 

iaa! 

is replaced by 

J dE b' ia {E) b ia (E) ^jj(E, a) if/(E, a) + c.c.} (35) 

ia 

and similarly for an analogous term in (ip\T\ip). Now, one has 

J2K a (E)b ia (E) = l -d E Y,bUEf = (36) 

i i 

by fl33|) . Therefore the only bi a dependent term in AT 2 is 

E / dE\b' ia {E)\^{E,a)\ 2 . (37) 

ia 

This becomes minimal if and only if b' ia = 0, i.e. 

bia(E) = Ci a , (38) 

with the real constant q q satisfying, by fl33|) . 

E(Q a ) 2 = 1 , a = ± . (39) 

i 

Inserting this into ( )32l) gives finally 

n = ^y d J E^ , {|^,+)(^ , ,+l + l^,-),^',-!} . (40) 

From ([3]) and ([1]) one now obtains the corresponding operators lit and T. Thus, in this 
one-dimensional example, covariance under time reversal and reflection plus minimality 
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of AT lead to a unique clock time operator and a unique associated temporal probability 
distribution. 

Interestingly, when the potential V in ( 1301) allows the application of scattering 
theory, n in (HOj) can be re-expressed in terms of the M0ller operators £l± and complex 
reflection and transmission coefficients, T(k) and R(k). We denote by \k) momentum 
eigenstates. Reflection invariance gives oVL±o = £l± and R(—k) = R(k), T(—k) = T(k). 
The interpolation property HQ± = Cl±H shows that the two vectors 

(ti + + n_)(\k) + \-k))/V2 

i(Cl + + Cl-)(\k)-\-k))/V2 (41) 

are eigenvectors of H which are symmetric and antisymmetric, respectively, and 
invariant under <j0 since 0f2±G = Cl T . Hence the vectors in (|4ip are real multiples 
of \E, ±). To determine the respective multiple we use the identities 

ti + + n_ = n_(s + i) = n+(i + &), (42) 



S\k) — T(k)\k) + R(k) \ — k) . (43) 



Thus 



S(\k)± | - fc» = (T(k) ± R(k))(\k) ± | - fc» . (44) 

This implies that here T(k) ± R(k) are pure phases. With fj42|) the two vectors in (j4ip 
become 

(1 + T(kj + R(k)) Q+(\k) + I - k))/V2, 

i(l + T{k)+R{k))Cl + {\k)-\-k))/V2. (45) 

Hence, from normalization one finds 

i(i=Fi)/2 ./m/Wi, A A _ 

|fl - ±) = |i + r( t) + fi(t)| ( tu + "- ) " t)± '-*» / ^- (46) 

Using ( 144")) this can be re-expressed by T, R, and f2+ as 



^■^"^fti^Hi^ 1 *^ 1 -*^- (47) 

Inserting this into ( )40l) and introducing as integration variable, J5/7i = ft k 2 /2m, one 
obtains an expression for Ho with coherences between \k) and | — A;). For a free particle 
in one dimension (T — 1, R — 0, f2 + = 1) one obtains 

n^ ) = / dk / dk'VWilfyik'] + | - fc}<-A/|} (48) 

so that in the free case there are no coherences between positive and negative momenta. 

Example 2: We consider the three-dimensional case of the Hamiltonian 

H = £ + V(r) (49) 
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and assume a purely continuous spectrum. H is invariant under rotations, reflection 
a i at the X2 — £3 plane, i.e. <fi — > 71 — <ft in spherical coordinates, and under time 
reversal. In the following we assume only invariance of T and its probability distribution 
operator under rotations and covariance under the combined action a\®\ this means the 
invariance of flo under these operations. Again this implies the existence of a unique 
time operator and its probability distribution, as shown now. 

The eigenfunctions of H can be chosen in the form J 'Ei(r)Yi m ($ <f>) , f Ei real. We 
denote these states by \Elm), with the normalization 

(Elm\E'l'm'} = 6(E - E')5 iv 5 mm> . (50) 
The action of <7iG is given by 

a^Elm) = (-l) m \Elm) . (51) 
Invariance under rotations and under <3"iB implies that n can be written in the form 

ft o = 2^Z)/ dEdE'J^bimiiE) \Elm){E'lm\ b ilm (E'), (52) 

i Ira 

with real functions bu m . In principle, the allowed set of / and m in the sum might depend 
on E, E'. However, if scattering theory can be applied this is not so, as seen further 
below. Note the diagonality in l,rn. The condition for the functions bu m , arising from 
the normalization of the probability, now reads 

(bu m (E)) 2 = 1 for each I, m . (53) 

i 

Equations ( l34l) - (|371) remain the same, only with a replaced by Im, and therefore 
minimality of AT implies 

flo = — VV I dE dE'\Elm){E'lm\, (54) 
2nn ' ^-^ J 

i Ira 

where (x.\Elm) = /Bz(r , )F; m ('i9, 0), with f Ei real. Thus also in this example the time 
operator is unique and there are no coherences between different lm values. 

Using scattering theory, (1541) can again be expressed by the M0ller operators Vt± 
and partial waves shifts, 8\. We write \Elm)o when V = 0. Then 



(x\Elm) = ^J^j l (kr)Y lrn (ti,<S>) (55) 

and, from rotation symmetry, 

S\Elm) = e 2iSl \Elm) . (56) 

Cl± commutes with rotations and a, and again Q&±Q = £l T . From H Cl± = &±Hq 
it follows that (Cl+ + (l-)\Elm)o is an eigenstate of H,L 2 ,L^, and also of &iQ with 
eigenvalue (— l) m . Therefore it is a real multiple of \Elm). To calculate this multiple 
we use (J4"2l) and obtain 

(0+ + n_)\Elm) = (1 + e T2iSl ) Vl ± \Elm) Q . (57) 
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Taking scalar products one obtains, from the unitarity of Cl 



\ Elm ) = 1 1 , 1 2ig ■ (0+ + 6_)|£7Zm)„ (58) 
1 + e zl0i 



and, from ( 1571) . 



I + e T2i5! 

II + e 2ld >\ 



\ Elm ) = H a&i n±\Elm)o • (59) 



Inserting this into fl54l) for il , using /c as integration variable, where E/K — hk 2 /2m, 
and the states 

\klm)o = H\/kJm\Elm)o, (60) 
which are normalized to 5(k — k')8u>8 mm >, then 

t- POO /"OO 

n = — V / dk / dk' 

1 + e '2i8,(k) 1 + e 2i5 ! (fc') 

VW |i + e a*,(fc)| |1 + e ^(fc0| l fe Moo(^m| fi 1 ,. (61) 
The free particle case results if one puts f2 + = 1 and 5j = 0. 

4. Application to arrival times 

Kijowski [13] considered, in three dimensions, a free particle of mass m coming in from 
X\ = — oo and determined the distribution of times of arrival at a plane perpendicular 
to the x\ axis, e. g. at X\ = 0, by using physically motivated postulates. These were: 
covariance, normalization, invariance under all Galilei transformations which leave the 
plane x\ — invariant, invariance under the combined action of time reversal and space 
reflection a: x — > — x, and minimal variance. His result, denoted here by fl t y , for the 
probability distribution of particles with only positive momenta in the X\ direction and 
for arrivals at the plane X\ — 0, can be written as 

^ poo /*oo /*oo 

^(^-^h/s- fc 2? A; 3 >(^, fc 2 , fc 3 | • (62) 

In the following, the methods of the previous sections will be used to generalize this 
result in two ways. First, invariance under Galilei transformations can be replaced by a 
weaker condition, and second, the restriction to particles coming in from one side can be 
lifted. We assume covariance, invariance under translations which leave the plane X\ — 
invariant, invariance under the space reflection a: x — > — x, time reversal covariance, 
and minimal variance. These symmetry assumptions are taken over from the classical 
case. Indeed, if one has an ensemble of free classical particles and reflects each particle 
trajectory by a, then the resulting reflected ensemble has the same distributions for 
arrivals at the plane X\ = as the original ensemble. In addition, one can go over 
to the time reversed trajectories, defined by Xe{t) = x(—t). Then, if a particle arrives 
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at the plane at time t', the corresponding particle on the time reversed trajectory will 
arrive there at time —if . This implies that for the time reversed ensemble the arrival 
time distribution is obtained from the original one by replacing t by —t. Note that the 
arrival location, i.e. the plane X\ = 0, is singled out here by the required reflection 
invar iance under a. 

Under the above assumptions it will now be shown that ITj t , the operator for the 
probability distribution of times of arrivals at the plane xi, is given by 

n? o = 0(«i) N 1/2 ^i)N 1/2 0(£i) 

+ 9(-v 1 )\v l \ 1/2 s(x 1 )\v l \ 1 / 2 9(-v 1 ), 

f[f t = exp{iH t/h}flf exp{-iH t/h}, (63) 

where V\ and X\ are components of the velocity and position operator and 9{x) is the 
usual step function (9(x) = for x < and 9(x) = 1 for x > 0). Note that classically 
the first term on the right-hand side corresponds to 9{v\) v\ S(xi) which classically gives 
the number of particles crossing the plane x\ per second from the left. Similarly for the 
second term, and their sum corresponds to \v\ \ S(xi). Thus in this case the arrival time 
distribution operator can be regarded as a particular - positive - quantization of the 
classical current. The usual quantum mechanical probability current is not necessarily 
positive, not even for states with only positive momenta in the x\ direction |15j . 
Inserting in (1631) two complete sets of momentum eigenstates and using 

(k^xOlk') = ±5(k 2 - k' 2 )5(k 3 - k' 3 ), (64) 

fl^ t can be written as 

fc roc roc roc 

flf = _ / dk 2 dk 3 dk X dk ljm-y?)tl2m 

1 27miJ_ oc J J 

\/k 1 k' l {\ki,k 2 ,k z )(k' 1 k 2 h\ + | - k x , k 2 , k 3 )(-k[, k 2 , k 3 \}. (65) 

For states with ki > this reduces to fl62l) . Note the absence of coherences between 
positive and negative k\ values in the above expression. For dimensions d other than 
three an analogous result holds. For d = 1 it becomes, under a change t — > —t, identical 
to the clock time expression in (f48l) . 

To show (1651) we introduce the notation 

|k, +) = (|fci, k 2 , k 3 ) + | - fci, k 2 , k 3 ))/V2, 

|k, -) = i {\k x , k 2 , k 3 )-\- h, k 2 , k 3 ))/V2, (66) 

which is a 5- normalized basis and invariant under <jO since <rO|k) = |k). The latter 
also implies that (kllljQlk') must be real. One now has 

(k,+|n^ |k',+) = (k,-|n^ |k',-) = (k|fr; |k'>, 

(k,+|fl^ |k',-) = 0. (67) 
Translation invariance in the x\ = plane gives the general form 

(k, +\tlf \W, +) = 5(k 2 - k' 2 )5(h - K)M(h, k 2 , k 3 , k[), (68) 
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where, for fixed k 2 and k 3 , M(-, k 2 , k 3 , •) is a positive definite kernel and one can write, 
with real functions bi(ki, k 2 k 3 ), 

M(h, k 2 , k 3 , k' 1 ) = 7 ^-Y] kih; k 2 , hMkt, k 2 , k 3 ) (69) 



(cf. Appendix A. 3 for more mathematical details). From this and (1671) we obtain, with 
|k',±> = |£^ 2 ,fc 3 ,±>, 

poo poo poo poo 

^(A: i; A: 2 ,fc 3 ){|k, +)(k',+| + \k,-)(k',-\}b i (k[;k 2 ,h). (70) 

Note again the diagonality, this time in k 2 and k 3 . Since Ily commutes with translations 
and hence with the momentum operators P 2 and P 3 , the normalization condition for 
the probability distribution 

— (il)\e i6t / h ILi e- iAt / n \ip) (71) 

for states with (k|^) vanishing at k\ = 0, gives the condition 

h ^ r „ „ f°° f h 



/poo 
dkidk 3 J dkidk[ 5 



x fc 2 , £3) &i(&i; &2, ^3) j^+^i^s) vp+iKhh) + (+ ^ - 

= X) / ^ 2 rfA;3^A; 1 -i6 J (A; 1 ;A; 2 ,fc3) 2 {|^ + | 2 + |^-r} = l, (72) 
where ^>±(k) = (k, ±\ip). Since \ip) can range through a dense set this yields 



^ Mfci^,fc 3 ) j =ij fci>Q _ (73) 

For states |^) with ^±(k) vanishing at k\ = one easily calculates that (ip\T\ip) is 
independent of hi and that for AT 2 the only hi dependent term is of the form 



2 



J dk 2 J dk 3 1^ dh^Y, ( d ^) 2 i\^\ 2 + \^-\ 2 } 



(74) 



where again the reality of bi is crucial. For AT 2 to be minimal this implies di(bi/ \fk\) = 
or 

—=bi(h] k 2 , k 3 ) = d(k 2 , k 3 ) . (75) 
Vh 

From ( I73|) it follows that 

J2( Cl (k 2 ,k 3 )f = l. (76) 

i 

Inserting ( I75|) into (!70|) and using ( IT6]) then gives 

^ /»oo /»oo /»oo /"CO 

flf n = / dk 2 / d/c 3 / dfcj / dk\ 

1 2ixmJ_ 00 J_ 00 J J 

V^llk, +){k[k 2 k 3 , +| + |k, -){k[k 2 k 3 , -|} . (77) 
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Inserting for |k, ±) from (15^]) one obtains the expression in f)65p for t = 0. The operator 
fl^ t is then obtained from (Q. Since Tlf commutes with P 2 and P3, only P x 2 remains 
and this yields ( 1651 ). 

The absence of coherences between positive and negative k% components in (1651) 
is due to the assumed reflection invariance. Recall that invariance under Galilei 
transformations as in [13] was not used here, only invariance under translations of the 
plane X\ = is assumed. 

If instead of the plane X\ = one considers arrivals at the plane X\ = a, the 
corresponding probability distribution operator is obtained by a spatial translation, 
which leads to an additional factor exp {i (k± — k[)a} in ( 1651 ). 



5. Interactions: conditional and operator-normalized arrival times 

In the derivation of the general form of covariant time operators the normalization to 1 
of the probability distribution is an essential condition. However, a free particle in two 
dimensions will arrive at a finite interval with a probability less than 1, and similarly 
in three dimensions. For a particle in a potential, also the arrival probability at a plane 
may be less than 1, due to scattering, and moreover, there may be several arrivals. 
Experimentally, in such a case one may measure the first time of arrival at a given 
position for a large number of replicas of the system and determine their distribution 
with respect to the total number of replicas. This distribution is, in general, not 
normalized to 1 since there may be systems with no arrival occurring. Dividing by 
the total probability for a first arrival or, equivalently, considering only those systems 
for which an arrival has actually been found, one obtains the conditional first-arrival 
probability distribution. Although this is, by construction, normalized to 1, it is not 
bilinear in the state vector and not the expectation of an operator. Alternatively, one 
may consider operator normalization as in [T6l \T7\ [TB] and then apply the techniques of 
Section [2j The physical relevance and interpretation of operator-normalized distribution 
in terms of a modification of the initial state has been pointed out in [IB] . 

The operator for the (possibly non-normalized) distribution of first arrivals is again 
denoted by fl/ 4 and the operator for the total probability for arrival at a given location 
by N. One has 

/oo 
dt tif . (78) 
-00 

For a given state \ip) the conditional distribution, denoted by II ^ for the times of first 
arrivals is the given by 

n a? = my> (79) 

(W\i>) 

which is evidently not bilinear in \ip). Since N is a positive operator, N~ 1 ^ 2 exists on 
states on which N does not vanish and which can be taken care of by a projector. If 
the external potential remains finite or is not too singular one expects from tunneling 
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that for any state there is a finite arrival probability, i.e. that N\ip) will never vanish. 
One can therefore define the operator fl^ t by 

Tlfi t = N~ l/2 fl t A N-V 2 . (80) 



Its time integral is then the unit operator and therefore the distributions determined by 
11$ t are normalized to 1. If iV and 11$ 4 are known one obtains IT/ 1 from (|80|) . 

To the normalized operator distribution fl^ t one can apply the procedures of 
Section [2] if additional symmetries hold. Physically motivated symmetries of fl t A carry 
over to N and thus also over to tl^t- More tricky is the requirement of minimal variance 
under given symmetry conditions. Physically, one would require minimal variance for 
the conditional distribution, i.e. that fl t A is chosen in such a way that n has minimal 
variance for each among all candidates satisfying the symmetry requirements. From 
f J79|) one has 



'IN 1 / 2 



■0 



IN 1 /* 



(81) 



If one defines the normalized vector \iPn) by 

Vjv = -=l ^- (82) 

\\N^)\\ 

one can write 

ILif = (MN-^ilfN- 1 ' 2 ^) = (^in^lVAr). (83) 

For given N, the correspondence between and \iPn) is one to one and therefore, if 
IT/ 1 leads to minimal variance the same is true for fl$ 4 , as long as N is kept fixed. 
Therefore, one may look for a normalized operator distribution satisfying required 
symmetry conditions and minimal variance, as in Section [2j If in addition the operator 
N for the total probability for arrival at a given location is known, one then obtains the 
physically interesting conditional probability distribution from Eqs. (|80|) and ( I79|) . 

To exemplify this, we consider one-dimensional motion in a bounded potential. 
The latter is assumed to be invariant under space reflection about the origin, and we 
consider an ensemble of particles which originally came in from ±oo. For an ensemble 
of classical particles there are two types of individual trajectories. A particle coming 
in from -oo (+oo) will either be reflected before it reaches the origin or, if it reaches 
the origin, will then continue to +oo {-oo), by the symmetry of the potential. We 
consider first arrivals at the origin. The corresponding classical arrival time distribution 
has certain symmetries. Indeed, from the reflection invariance of the potential, in the 
ensemble obtained by space reflecting each trajectory of the original ensemble has the 
same distribution for arrivals at the origin as the original ensemble. Going over to the 
time reversed trajectories one concludes, as in Section H] for the time reversed ensemble, 
that the arrival time distribution is obtained from the original one by replacing t by — t. 

These symmetry properties are now carried over to the quantum case by demanding 
corresponding symmetries for II f, i.e. invariance under space reflection and covariance 



Symmetries and time operators 



16 



under time reversal. This implies the corresponding properties for the operator- 
normalized tl^ t and in particular invariance of 11$ under space reflection and time 
reversal. This then is the same situation as in Example 1 of Section [3] where a clock 
operator with the similar symmetry properties was constructed. Therefore, for n$ on 
obtains the same result as in (HDD, where, as a reminder, \E, ±) denote eigenstates of H 
such that (x\E, ±) are real wavefunctions which are either symmetric or antisymmetric. 
/,From this expression and Eqs. (j7|) - ffTTl) one then obtains fi^. 

To determine the corresponding conditional probability distribution N is needed, 
the operator for the total arrival probability at the origin. An at least approximate 
expression for N may be obtained by modeling the detection process. The simplest way 
to do this is by a complex potential at the origin or by a fluorescence model (cf. [19] 
for a review of these models). From the N calculated in this way and from the above 
result for fl^t one then obtains an expression for the conditional distribution of arrivals 
at the origin. The complex potential and fluorescence models also allows the calculation 
of an approximate expression for the conditional arrival distribution, and one may then 
compare the two expression, in particular whether the abstractly obtained operator- 
normalized distribution together with the approximate result for N leads to smaller 
variances than the distributions derived from the models. 

This example can be carried over to three dimensions and a Hamiltonian of the 

form 

H = P 2 /2m + V(\X 1 \) , (84) 



using techniques similar to those described in Appendix Appendix A. 3 



6. Discussion 



We have derived the most general form of covariant time operators with a normalized 
probability distribution. The result has been applied to clock time operators and to time- 
of-arrival operators. To arrive at a unique operator, the multitude of possible operators 
has been restricted by physically motivated symmetry requirements, in particular time 
reversal invariance and, to obtain an optimally sharp observable, by the requirement of 
minimal variance. 

Our general result has allowed us not only to weaken the assumptions of Kijowski 
[T3] but at the same time generalize his arrival-time distribution for free particles arriving 
at a plane. It had been required in [13] that the distribution should be invariant under 
Galilei transformations which transform the plane into itself. As shown here, Galilei 
transformations are not needed and translation invariance is sufficient. Moreover, the 
restriction that the particles should come in from one side of the plane only has been 
lifted. 

A free particle in three dimensions will arrive at a location of finite extent only with 
a probability less than 1, and also for a particle in a potential the arrival probability at a 
plane may be less than 1, due to scattering, and moreover, there may be several arrivals. 
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In an actual and repeated repeated measurement one can consider only those outcomes 
for which at least one arrival has actually been found at some time and then in addition 
consider only the distribution of the times of the first arrival. This amounts to a division 
by the total arrival probability and determines the (normalized) conditional first-arrival 
probability distribution. However, the latter is not bilinear in the state vector, due to 
the division. As an alternative, we have applied operator normalization to the original 
distribution, and then one can apply the general results. The distribution thus obtained 
leads to the conditional distribution if the total arrival probability is known. 

It should be pointed out that the probability distribution associated with an 
observable is in general of more direct experimental relevance than the operator of 
the observable itself because in most circumstances first the distribution of individual 
measurement results for a given state is obtained and then the expectation value is 
derived from that. For a particle in a potential, the total arrival probability at a plane 
may be less than one. In such a case there may still be a non-normalized temporal 
probability distribution operator but not necessarily an associated time operator. As 
has been pointed out above, the conditional probability distribution is in general not 
bilinear in the state vector and so cannot be used directly to construct a time operator. 

In order to be meaningful, any experimental procedure purported to be the 
realization of a measurement of arrival times or of a quantum clock has to satisfy 
certain requirements, and the so constructed ideal observables provide guide-lines for 
these requirements. In this context, an analysis of this kind for the ideal arrival time- 
of-arrival distribution of [13] has been carried out in terms of an operational quantum- 
optical realization with cold atoms (cf. [19] for a review). 
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Appendix A. Mathematical details 

Appendix A.l. Nondegenerate energy spectrum 

In this appendix, mathematical details are given to justify some more formal arguments 
in the main text. To illustrate the crucial points we first treat the case of a Hamiltonian 
H with a simple absolutely continuous spectrum from [E ,oo) or (—00,00), with 
no degeneracy. The generalized eigenvectors can then be denoted by \E), with 
normalization (E\E') = 5(E — E'). In the Hilbert space, "H, of the Hamiltonian we 
denote by S(H) and T>(H) the subspaces of vectors € H such that (E\ip) = ip(E) 
is in the Schwartz space S(R) or T>(R), respectively. Then n = —i[H,F ]/h from 
01]) is a positive semi-definite continuous sesquilinear form on S(H) or T>(H) and, by 
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continuation, on S(R) or D(R). Hence, by the kernel theorem, it defines a distribution 
on 5([R 2 ) or V(R 2 ), respectively, formally denoted by 

U (E,E') = (E\fl \E') . (A.l) 
Let ||<?i) G T>{H) , i — 1, ■ ■ • | be maximal set (finite or infinite) such that 

(</i|n |&) = 6 tj . (A.2) 
We put gi(E) = (E\gi) and define h G S'(R) by 

bi(E) = J dE'U (E,E') gi(E') . (A.3) 

Then 

(g a \bi) = J dEg~jE)b i {E) = S ai (A.4) 
and, therefore, 

Y^(9a\bi) (bilge) = 2nh <W ■ (A.5) 

i 

With £{gi} the linear span of the gi(E) 's one therefore has, on £{gi} x 

(E\U \E') = -Lj^^ipJJ . (A.6) 

j 

If A/" C «S(R) or C P(R) defines the null space of (the distribution defined by) n , then 

£{gi} U Af span S(R) or P(R), respectively, and therefore, by continuity, ( 1A.6I) holds 
also on <S(#) x S(H) and X>(if) x V(H), respectively 

Exploiting normalization. From Eqs. ([SD and (I A . 6 [) one obtains, for a normalized 



1= / A^le-^IIoe^l 

,/ Z7T/1 



1 



27T 



dEe lt,T bi(E)i/}(E) 



(A.7) 



This implies that, for each z, the inner integral, considered as a function of r, is square 
integrable, and so is its Fourier transform. The latter, as a distribution, is given by 
bi(E)ip(E) and hence is equivalent to a square integrable (i.e. L 2 ) function. Hence 
bi(E) is not only a distribution but even locally an L 2 -function. From (1A.7J) one then 
has by Parseval's theorem 

J2 fdE \h(E)i)(E)\ 2 = 1 . (A.8) 

i J 

Since this is true for any normalized \ip) G S(H), this implies 

H E ) h{E) = 1 for almost all E. (A.9) 
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Conversely, if a sequence of locally L 2 -functions satisfies (IA.9j) . then II defined by f]A.6|) 
gives rise to a covariant time operator through Eqs. ([5]) and ([I]). 

Exploiting the second moment. We assume that for a dense set of \ip) in S(H) the 
second moment exists, i.e. 

/ ^-t 2 (V|e-^ t/R noe^ t/B |^) < oo . (A.IO) 
J 2nn 

With (IA.6I) this can be written as 

H 2 J2 [ drr2 / dEe iET b~(E)tP(E) . (All) 

This implies that, for each i, the inner integral considered as a function of r is not 
only square integrable but is also in the domain of the operator "multiplication by 
t". Therefore, by standard results, the Fourier transform of this function, given by 
bi(E)ij)(E), is not only in L 2 but is also absolutely continuous and has a derivative 
which is also in L 2 . Again by Parseval's theorem one has 

second moment = ti 2 ^ f dE d E (bi(E)ij;(E)\ < oo . (A.12) 

i J 

Since is from the dense set S(H), h must also be absolutely continuous and its 
derivative exists almost everywhere and is locally L 2 . Using (1A.9j) one obtains 

second moment/^ 2 = (A. 13) 

JdE\d E i)\ 2 + fdEj2\K^\ 2 + 2Re JdE^bTbi^'il; (A.14) 

In a similar way one finds 

(iP\f |V) = JdE4>ihd E ip + J dE\iP\ 2 J2 b i ihd Ek ■ (A.15) 

We note that if the hi 's are real then (1A.9j) implies bi(E) b' { (E) = for almost all E 
and thus for real hi 's the last term in ( 1A.14I) as well as in ( 1A.15I) vanishes. 

Appendix A. 2. Extension to degenerate energy spectrum. 

In case of a discrete degeneracy label we use the normalization as in (flBl where for 
simplicity we assume the same degeneracy for each E. (The general case can be 
treated by direct integrals of Hilbert spaces). Then T1(E, E'), gi(E) and bi(E) from the 
nondegenerate case are replaced by U(E,a, E',a'), gi(E,a) and bi(E,a), respectively, 
and ( 1A.6j) can be written as 



(E, a\fl \E', a') = ^— y>(£, a) h{E', a') . (A.16) 

i 

By the same technique as above the normalization condition now implies 

J2bi(E,a)b i (E,a') = 5 aa , . (A.17) 
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Finiteness of the second moment implies that bi(E, a) is absolutely continuous in E for 
each a and its derivative is locally L 2 . 

In case of a continuous degeneracy label we use the normalization 

(E, X\E'X') = 5(E - E') S(X - A') (A.18) 

This can be reduced to the discrete case by choosing a complete orthonormal set 
{$ a (A), a = 1, 2, • • ■} and putting 

\E, a) = J dX $ a (A) \E, A) , a = 1, 2, • • • (A.19) 

It can also be treated directly by replacing the functions gi{E) by g%{E, A) and bi(E) by 
bi(E, A). Then ( 1A.17I) is replaced by the equation 

HE, A) bi(E,X') = 5(X - A') (A.20) 

i 

valid for almost all E. 

Appendix A. 3. The case of partial translation invariance 

To illustrate some mathematical details for Section H] we consider for notational 
simplicity the two-dimensional case. We search for a covariant probability distribution 
operator, denoted by flf, which is invariant under reflection a : x — > — x, under 
translations which leave the line x\ — invariant, as well as under time reversal 0, 
and which has minimal variance. Let H\ be the one dimensional Hamiltonian 

H x = Pi /2m (A.21) 

in the Hilbert space "Hi = L 2 (R) and denote by \Ei,±) the real symmetric resp. 
antisymmetric eigenstates of Hi, similarly as in Example 1 of Section [31 The eigenstates 
of H are then 

\Ei,±)\k 2 ) . (A.22) 

The assumed invariances imply the corresponding invariances of fl^ and of F^. The 
expectation value of the latter gives the probability of finding a measured time value in 
(— oo, 0), and so < Fq < 1. If \ip) E H, the formal expression 

mh)} = (k 2 \iP) (A.23) 

can be regarded, for each fc 2 , as a vector in "Hi, i.e. as a vector function of k 2 with values 
in %, and for j^i), l^) G % one has 

dk 2 ^i{k 2 )\i} 2 {k 2 )) Hl = (Vi|-0 2 > (A.24) 

where the index indicates the scalar product in %i. Each bounded operator A which 
commutes with P 2 satisfies 

(k 2 \A\k' 2 ) = S(k 2 -k' 2 )A(k 2 ) , (A.25) 
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where A(k 2 ) acts in Hi, and thus 

(k 2 \A\i>) =A(k 2 M(k 2 )) . (A.26) 
In a rigorous way this can be expressed in terms of direct integrals of Hilbert spaces, 

n= [ dk 2 U{k 2 ) , where U{k 2 )=U x , (A.27) 
J® 

with |^) corresponding to a vector-valued function \ip(-)), with \ip(k 2 )) G H(k 2 ), and 
A corresponding to an operator- valued function A(-) satisfying flA.26j) . If A is positive 
and bounded by 1 then so is A(k 2 ), for almost all k 2 . 

We now take for A. Since Hi acts in each 7-L(k, 2 \ the action of fl^ is given by 



fl£(k 2 ) = -1 [H u F A {k 2 )\ (A.28) 

considered as a bilinear form in each 1-L{k 2 ). For almost all k 2 , IlQ(k 2 ) is a positive 
bilinear form satisfying invariance under a and B. Hence each 11^ (k 2 ) is given by fHUj) . 
(with E replaced by Ei). This finally implies 

n^ = ^y dE 1 dE[{\E 1 ,+)(E 1 ,+\ + I*?!, -|} (A.29) 

or, equivalently, 

= dk 2 dE 1 dE[ 

[^\E 1 +)\k 2 )(k 2 \(E' 1 ,+\^) + (+**->}. (A.30) 

In three dimensions, the procedure is analogous. Now ( 1771) results by a change of 
variable. 
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